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Abstract 

The  current  effort  develops  and  demonstrates  the 
application  of  high  resolution  turbulence  modeling  to 
flight  mechanics  and  aeroelasticity  of  air  vehicles  at  flight 
conditions  where  the  vehicle  is  experiencing  massively 
separated  flow  fields.  The  effort  has  both  a basic 
research  component  to  aid  in  developing  the  method  and 
an  applied  component  where  the  method  is  used  to 
demonstrate  an  ability  to  simulate  current  DoD  aircraft 
issues  in  flight  mechanics  and  aeroelasticity.  The  high 
resolution  turbulence  method  is  a hybrid  Reynolds 
Averaged  Navier-Stokes  (RANS)-Large  Eddy -Simulation 
(LES)  method  introduced  by  Spalart  et  al.  in  1997  called 
Detached-Eddy  Simulation  (DES)  implemented  in  an 
unstructured  Navier-Stokes  solver.  Cobalt. 

In  the  basic  research  component,  DES  has  been 
applied  to  an  Aerospatiale-A  airfoil  at  an  angle  of  attack 
of  13.3  degrees  and  a Reynolds  number  of  2 million.  The 
project  is  called  DESFOIL  and  simulates  laminar-to- 
turbulent  transition,  adverse  pressure  gradients, 
streamline  curvature,  and  boundary  layer  separation  of  a 
3-D  airfoil  strip.  This  study  is  in  the  early  stages  of 
developing  a baseline  for  RANS  and  DES  computations. 

DES  has  also  been  applied  to  flight  mechanic  and 
aeroelasticity  problems  of  DoD  air  vehicles  to 
demonstrate  the  utility  of  DES  and  also  discover  some  of 
the  nonlinear  mechanisms  causing  these  flight  issues.  The 
applications  studied  include  the  F/A-18E  forced  motion 
about  the  roll  axis  and  one  degree  of  freedom  simulation 
of  abrupt  wing  stall  (AWS),  the  F/A-18C  at  conditions  of 


Boeing 
Commercial 
Airplanes,  Seattle, 

I WA 

Philippe.r.  Spalart 
@boeing.com 

i 
i 
I 

i 

tail  buffet,  and  the  ARGUS  missile  at  conditions  where  it 
experiences  coning  motion.  I 

i 
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1.  Introduction  I 

i 

f 

This  work  focuses  on  multidisciplinary  applications 
of  Detached-Eddy  Simulation  (DES),  principally  flight 
mechanics  and  aeroelasticity.  Specifically,  the  lateral 
instability  (known  as  abrupt  wing  stall)  of  the  pre- 
production  F/A-18E  is  reproduced  usingjDES,  including 
the  unsteady  shock  motion.  A single  degree-of-ffeedom 
calculation  is  performed  as  well  to  demonstrate  the  onset 
of  the  wing  drop.  DES  is  applied  to  the  F/A-18C  at  a 
moderate  angle  of  attack  to  reproduce  the  vortex 
breakdown  leading  to  vertical  stabilizer  buffet.  Unsteady 
tail  loads  are  compared  to  flight  test  data.  | 

Previous  DoD  Challenge  Project  work  has 
demonstrated  the  unique  ability  of  the  DES  turbulence 
treatment  to  accurately  and  efficiently  predict  flows  with 
massive  separation  at  flight  Reynolds  numbers.  DES 
predictions  are  obtained  on  unstructured1  grids  using  the 
Cobalt  code,  an  approach  that  can  accommodate  complete 
configurations  with  very  few  compromises.  A broad 
range  of  flows  has  been  examined  in  previous  DoD 
Challenge  Project  work,  including  aircraft  forebodies, 
airfoil  sections,  missile  afterbody,  vortex  jbreakdown  on  a 
delta  wing,  and  the  F-16  and  F-15E  at!  high  angles-of- 
attack.  All  DES  predictions  exhibited  a moderate  to 
significant  improvement  over  results  | obtained  using 
traditional  Reynolds-averaged  models  ana  often  excellent 
agreement  with  experimental/flight-test  data  is  observed. 
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DES  combines  the  efficiency  of  a Reynolds-averaged 
turbulence  model  near  the  wall  with  the  fidelity  of  Large- 
Eddy  Simulation  (LES)  in  separated  regions.  The 
development  and  demonstration  of  improved  methods  for 
the  prediction  of  flight  mechanics  and  aeroelasticity  in 
this  DoD  Challenge  Project  is  expected  to  reduce  the 
acquisition  cost  of  future  military  aircraft. 

The  F/A-18E  flight  mechanic  simulations  are  in  the 
second  full  year  of  development.  The  previous  year, 
static  calculations  were  made  of  a full  and  half  span 
model  at  conditions  where  it  experiences  AWS.  These 
simulations  compared  very  favorably  with  experimental 
data.  This  year’s  effort  incorporated  dynamic  motion  of 
the  vehicle  along  the  roll  axis  to  mimic  wind  tunnel  free- 
to-roll  experiments.  Two  different  pitch  angles  were 
examined  and  unsteady  data  was  obtained  and  compared 
to  the  experimentally  obtained  frequency  data.  The 
current  calculations  have  so  far  qualitatively  captured  the 
experimental  data.  This  application  has  made  great 
strides  in  demonstrating  the  utility  of  using  an 
unstructured  solver  and  DES  to  compute  the  critical 
nonlinear  aerodynamics  necessary  to  estimate  static  and 
dynamic  control  derivatives  of  fighter  aircraft. 

The  F/A-18C  tail  buffet  calculations  are  also  in  the 
second  year  of  development.  Simulations  were  performed 
on  a configuration  similar  to  the  F-18  High  Alpha 
Research  Vehicle  and  compared  to  flight  test  data. 

2.  Numerical  Method 

Solutions  were  computed  with  the  commercial 
version  of  Cobalt  developed  by  Cobalt  Solutions.1'1 
Cobalt  solves  the  unsteady,  three-dimensional, 
compressible  Navier-Stokes  equations  on  a hybrid 
unstructured  grid.  The  code  has  several  choices  of 
turbulence  models,  including  Spalart  Almaras  (SA)  and 
Menter’s  Shear  Stress  Transport  (SST)  Reynolds- 
averaged  Navier-Stokes  (RANS),  as  well  as  DES  versions 
of  SA  and  SST.  All  simulations  were  computed  on 
unstructured  meshes  with  prisms  in  the  boundary  layer 
and  tetrahedra  elsewhere  on  half-span  surface  geometries. 
The  computational  meshes  were  generated  with  the 
software  packages  GridTool121  and  VGRIDns.131 

For  simulation  of  turbulent  flows,  the  governing 
equations  are  suitably  averaged,  yielding  turbulent 
stresses  that  require  a model.  A Boussinesq 
approximation  is  invoked  in  the  momentum  equations  and 
the  turbulent  eddy  viscosity  (u,)  is  used  to  relate  the 
stresses  to  the  strain  rate.  The  turbulent  heat  flux  is  also 
modeled  using  a gradient-transport  hypothesis,  requiring 
specification  of  a turbulent  thermal  conductivity,  kt.  The 
Reynolds  analogy  is  applied  and  the  turbulent  heat  flux  is 
modeled  using  a constant  turbulent  Prandtl  number  of  0.9. 
Using  turbulent  eddy  viscosity  and  turbulent  conductivity. 


the  variable  n (u  + p,)  is  replaced  by  and  k is  replaced  by 
( k +k, ) in  the  governing  equations. 

DES  was  proposed  by  Spalart  et  al.14,51  The 
motivation  for  this  approach  was  to  combine  LES  with 
the  best  features  of  RANS  methods.  RANS  methods  have 
demonstrated  an  ability  to  predict  attached  flows  veiy 
well  with  a relatively  low  computational  cost.  LES 
methods  have  demonstrated  an  ability  to  compute 
separated  flowfields  accurately,  but  at  a tremendous  cost 
for  configurations  with  boundary  layers.  Spalart’ s DES 
method  is  a hybrid  of  LES  and  RANS,  which  combines 
the  strengths  of  both  methods. 

The  DES  model  was  originally  based  on  the  Spalart- 
Allmaras  one  equation  RANS  turbulence  model.  The 
wall  destruction  term  is  proportional  to  (v  / </)  , where  d 
is  the  distance  to  the  wall.  When  this  term  is  balanced 
with  the  production  term,  the  eddy  viscosity  becomes 
proportional  to  Sd2  where  S is  the  local  strain  rate.  The 
Smagorinski  LES  model  varies  its  sub-grid  scale  (SGS) 
turbulent  viscosity  with  the  local  strain  rate,  and  the  grid 
spacing:  vSgs  x SA2,  where  A = max(Ax,  A_y,  A z).  If  d is 
replaced  with  A in  the  wall  destruction  term,  the  S-A 
model  will  act  as  a Smagorinski  LES  model. 

To  exhibit  both  RANS  and  LES  behavior,  d in  the  SA 
model  is  replaced  by 

d = mm(d,CoesA) 

When  d« A,  the  model  acts  in  a RANS  mode  and 
when  d« A the  model  acts  in  a Smagorinski  LES  mode. 
Therefore,  the  model  switches  into  LES  mode  when  the 
grid  is  locally  refined. 

DES  was  implemented  in  an  unstructured  grid 
method  by  Forsythe  et  al.161  They  determined  the  CDES 
constant  should  be  0.65,  consistent  with  the  structured 
grid  implementation  of  Spalart  et  al.141  when  the  grid 
spacing  A was  taken  to  be  the  longest  distance  between 
the  cell  center  and  all  of  the  neighboring  cell  centers. 

A Newton  sub-iteration  method  is  used  in  the 
solution  of  the  system  of  equations  to  improve  time 
accuracy  of  the  point-implicit  method  and  approximate 
Jacobians.  In  the  calculations  presented  below,  a typical 
number  of  three  Newton  sub-iterations  is  used  for  all 
time-accurate  cases. 

2.1.  Summary  of  the  Proposed  Method. 

The  proposed  method  for  simulating  aircraft  at  flight 
Reynolds  numbers  in  conditions  of  massively  separated 
flow  is  as  follows: 

1.  Use  a time-accurate  unstructured-grid  solver 
with  moving  mesh  capability  to  allow  rapid  turn 
around  of  grids  on  complex  configurations  - the 
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solution  must  have  at  least  second-order  spatial 
and  temporal  accuracy. 

2.  Use  DES  as  the  underlying  turbulence  treatment 
to  obtain  accurate  unsteady  loads  and  mean 
quantities  - this  requires  a low  dissipation  solver. 

3.  Use  Adaptive  Mesh  Refinement  (AMR)  to 
improve  grid  resolution  in  critical  areas  with 
nonlinear  flowfield  phenomena. 

3.  Results 

Results  are  shown  for  both  the  basic  and  applied 
portions  of  the  challenge  project.  The  basic  project  is  a 
start-up  effort  to  investigate  laminar  to  turbulent  transition 
on  an  airfoil  shape.  In  the  applied  portion  of  the  study, 
two  full  aircraft  configurations  are  analyzed  for  flight 
mechanic  and  aeroelastic  phenomena  of  abrupt  wing  stall 
and  tail  buffet  and  then  the  two  configurations  are 
analyzed  as  a result  of  flight  test  support  requests. 

3.1.  DESFOIL. 

Prediction  of  complex  flows  that  include  laminar-to- 
turbulent  transition,  adverse  pressure  gradient,  streamline 
curvature,  and  boundary  layer  separation  remain  among 
the  most  challenging  for  turbulence  simulation  strategies. 
A prototypical  example  that  is  the  focus  of  the  present 
investigation  is  the  flow  over  an  airfoil  at  maximum  lift. 
Flow  regimes  are  sensitive  to  the  airfoil  geometry,  angle 
of  attack,  and  Reynolds  number  and  motivate  various 
hierarchies  of  simulation  strategies.  The  specific  flow  of 
interest  is  that  over  the  Aerospatiale-A  airfoil  at  an  angle- 
of-attack  of  13.3  degrees  and  Reynolds  number  of  2 x 
106,  corresponding  to  maximum  lift.  The  flow  has  been 
measured  in  separate  experiments  and  was  the  subject  of  a 
coordinated  set  of  investigations  through  the  LESFOIL 
project 

DES  and  Reynolds-averaged  Navier-Stokes 
predictions  have  been  obtained  of  the  flow  over  the 
airfoil,  with  the  objectives  to  date  being  to  establish  a 
baseline  upon  which  enhancements  to  the  predictive 
technique  can  be  assessed.  Shown  in  Figure  1 are 
contours  of  the  instantaneous  vorticity  in  four  planes 
along  the  airfoil.  At  x/C  = 0.4,  the  RANS  model  is 
retained  and  the  figure  illustrates  that  the  solution 
possesses  weak  spanwise  variation.  At  the  subsequent 
planes  a range  of  scales  is  resolved  as  the  flow  develops 
eddies  in  the  separating  shear  layer.  The  computations 
performed  to  date  have  successfully  demonstrated  that  the 
approach  of  handling  laminar-to-turbulent  transition  is 
numerically  feasible  and  relatively  accurate.  Further 
investigations  will  begin  the  process  of  incorporating 
eddy-seeding  strategies  into  the  simulations,  along  with 
substantial  grid  refinement  in  order  to  support  turbulent 


structures  within  the  boundary  layer  in  the  aft  region  of 
the  airfoil.  There  remain  very  significant  challenges  to 
the  modeling  strategy  that  will  require  substantial 
computational  resources,  in  turn  further  motivating  the 
need  for  HPC  resources. 

i 

3.2.  F/A-18E  Abrupt  Wing  Stall. 

During  envelope  expansion  flights  of  the  F/A-18E  in 
the  Engineering  and  Manufacturing  Development  phase, 
the  aircraft  encountered  uncommanded  j lateral  activity, 
which  was  labeled  “wing  drop”.  An  extensive  resolution 
process  was  undertaken  to  resolve  this  issue.  A 
production  solution  was  developed,  which  included 
revising  the  flight  control  laws  and  the  incorporation  of  a 
porous  wing  fold  fairing  to  eliminate  ( the  wing  drop 
tendencies  of  the  pre-production  F/A-18E/F.  The  wing 
drop  events  were  traced  to  an  abrupt  wing  stall  (AWS)  on 
one  side  of  the  wing  causing  a sudden  and  severe  roll-off 
in  the  direction  of  the  stalled  wing.  Development  of  a 
reliable  computational  tool  for  prediction  of  abrupt  wing 
stall  would  enable  designers  to  screen  configurations  prior 
to  building  the  first  prototype,  reducing  costs  and  limiting 
risks.  | 

The  F/A-18E  provides  an  excellent  testing  ground  for 
simulation  tools  due  to  the  large  amount]  of  experimental 
data  obtained.’7,8’  Previous  computational  research19’ 
focused  on  predicting  the  zero  sideslip  characteristics  of 
the  aircraft,  including  the  break  in  the  llift  curve  slope 
characteristic  of  AWS.  It  was  found  that  by  applying 
Detached-Eddy  Simulation  (DES)  to  this  problem  to 
predict  the  unsteady  shock  motion  seen  experimentally,  a 
better  mean  flow  prediction  could  be  obtained  compared 
to  industry  standard  Reynolds-averaged  (RANS) 
models.’10’ 

The  current  work  seeks  to  extend  the  past 
computational  successes  to  predicting  stability  derivatives 
(both  static  and  dynamic)  in  the  AWS  regime.  An 
unstructured  full  aircraft  grid  was  created  with  8.4*  106 
cells  by  using  a coarse  baseline  grid  and  then  using 
solution  based  mesh  adaptation  to  cluster  points  in  the 
separation  region  above  the  wing.  Both  Menter’s  SST 
RANS  model  and  Detached-Eddy  Simulation  were 
applied.  To  assess  the  accuracy  of  the  simulations, 
comparisons  are  made  against  experiments.  Normal  force 
vs.  angle-of-attack  is  plotted  in  Figure ' 2,  showing  the 
slope  break  in  the  experiments.  DES , shows  a better 
agreement  than  SST  RANS  in  this  case,'  as  was  seen  in 
previous  work. 

Calculations  were  also  performed  with  various  bank 
and  pitch  angles.  For  the  experiments  and  computations, 
the  pitch  angle  was  held  fixed,  and  the  model  rolled 
around  the  longitudinal  axis  of  the  aircraft.  This  leads  to 
a reduction  in  alpha,  and  an  increase  in  j beta.  Thus  the 
calculations  do  not  strictly  give  derivatives  with  respect  to 
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beta.  Figure  3 shows  rolling  moment  and  yawing  moment 
vs.  roll  angle  for  7°  pitch  angle.  The  agreement  to 
experiments  for  yawing  moment  is  quite  good,  since  this 
comes  mainly  from  the  vertical  stabilizer,  which  is  not 
separated.  The  agreement  for  rolling  moment  is  less 
accurate  since  the  unsteady  location  of  the  shock,  which 
separates  the  flow  on  the  wing,  is  challenging  to  predict. 
The  change  in  sign  in  rolling  moment  for  the  SST  RANS 
at  30°  bank  was  due  to  the  shock  on  the  down-turned 
wing  moving  forward,  decreasing  lift  on  that  wing. 

Calculations  have  also  been  performed  in  a forced 
oscillation  to  estimate  roll  damping.  A sample  flow 
visualization  from  a DES  simulation  illustrating  the 
separated  region  is  shown  in  Figure  4.  Also  shown  for  the 
DES  simulation  are  phase  averaged  plots  for  four  pitch 
angles,  each  derived  from  five  cycles  of  the  rolling 
moment  vs.  non-dimensionalized  roll  rate.  Stable 
behavior  (i.e.,  negative  slope)  is  seen  at  6°  and  7°.  At  8° 
there  are  some  strong  non-linearities,  while  the  9°  plot 
shows  regions  of  unstable  roll  damping. 

Clearly,  the  ability  to  computationally  predict  static 
and  dynamic  stability,  especially  at  transonic  and  higher 
Mach  numbers,  where  experimental  facilities  are  quite 
limited,  would  provide  a significant  increase  in  capability 
for  airplane  design  and  analysis.  Free-to-roll 
computations  are  currently  underway  to  compare  directly 
to  free-to-roll  experiments,  which  have  been  highly 
successful  in  correlating  to  flight  tests. 

3.3.  F/A-18C  Tail  Buffet. 

The  F/A-18C  simulations  were  conducted  to 
demonstrate  the  ability  of  the  method  to  reproduce  the 
aerodynamics  of  tail  buffet.  Tail  buffet  of  the  F/A-18C  is 
a fluid  structure  interaction  resulting  from  burst  leading- 
edge  extension  vortices  impacting  the  twin  vertical  tails 
and  was  observed  in  extensive  flight  tests  of  the  F-18 
FLARV.  At  realistic  flight  conditions  this  flow  field  is 
also  complicated  by  turbulent  flow  generated  in  the  post 
breakdown  region  surrounding  the  tails  and  in  the 
boundary  layer  of  the  vehicle.  Results  are  compared  to 
unsteady  tail  pressure  coefficient  data  and  vortex 
breakdown  locations  obtained  in  the  NASA  F-18  FLARV 
flight  tests.  Follow-on  studies  of  this  configuration  will 
incorporate  aeroelastic  tails  to  fully  simulate  the 
phenomena. 

All  F/A-18C  cases  were  run  at  30°  angle-of-attack,  a 
Mach  number  of  0.2755,  and  a standard  day  altitude  of 
20,000  feet.  The  resulting  Reynold’s  number  was  13 
million  based  on  the  mean  aerodynamic  chord  of  the 
aircraft  (12  ft).  The  baseline  grid  of  3.6  million  cells  was 
generated  with  VGRFDns.  Unsteady  SADES  turbulence 
model  simulations  were  performed  using  the  baseline 
grid.  A time-averaged  SADES  solution  was  used  to 
produce  an  AMR  grid  with  3.9  million  cells  by  following 


the  approach  outlined  in  Reference  4.  All  time-accurate 
simulations  were  run  for  over  10,000  iterations  with 
second-order  temporal  and  spatial  accuracy,  three  Newton 
sub-iterations,  and  a time  step  of  0.0005  seconds.  The 
chosen  time  step  results  in  a time  step  non- 
dimensionalized  by  the  freestream  velocity  and  mean 
aerodynamic  chord  of  0.0012.  This  characteristic  time 
step  was  found  adequate  in  previous  studies  of  vortex 
breakdown  and  massively  separated  flows.111’12’13,141 

Solutions  were  computed  using  the  SST,  SA,  and 
SADES  turbulence  models  to  determine  their  effect  on  the 
fiowfield.  Solutions  for  all  three  methods  were  computed 
using  the  same  grid,  time  step,  and  number  of  sub- 
iterations to  provide  a consistent  comparison.  Figure  5 a- 
c depicts  snapshots  of  solutions  for  each  method  with  the 
surface  colored  by  pressure  and  an  iso-surface  of  vorticity 
shown.  The  chosen  vorticity  level  for  the  isosurface  and 
the  pressure  colormap  are  held  fixed.  Although  the 
snapshots  are  not  necessarily  synchronized  in  time,  the 
overall  differences  are  striking.  The  SADES  solution 
(Figure  5c)  produces  a much  more  detailed  view  of  the 
simulation  since  it  is  able  to  capture  much  finer  fiowfield 
scales.  The  SST  (Figure  5a)  and  SA  (Figure  5b)  models 
are  unable  to  capture  the  proper  post-breakdown  behavior 
or  the  leading-edge  separation  regions  of  the  wing, 
horizontal,  and  vertical  tails.  It  is  also  apparent  that  the 
SST  LEX  vortex  pressure  footprint  on  the  surface  is 
significantly  different  than  either  the  SA  or  SADES 
solutions.  The  low  pressure  region  represented  by  a dark 
green  color  is  greatly  reduced  in  size  on  the  SST  solution. 
The  SADES  solution  is  also  capturing  the  vortical  sub- 
structures around  the  primary  vortex. 

A common  definition  of  vortex  breakdown  is  the 
location  where  the  streamwise  velocity  component  is  zero 
in  the  core.  The  coordinates  of  this  point  along  the  core 
were  tracked  in  time  for  each  of  the  methods,  SST,  SA, 
and  SADES.  Figure  6 depicts  the  time  histories  of  the 
three  methods  as  well  as  the  flight  test  and  experiment 
maximum  and  minimum  mean  values  of  vortex 
breakdown  presented  in  Reference  15.  Three  things  are 
obvious  from  Figure  6.  First,  the  amplitude  of  oscillation 
for  the  SST  and  SA  models  is  almost  negligible  compared 
to  the  SADES  simulation.  Second,  the  SST  solution 
predicts  breakdown  far  upstream  of  the  flight  test  or 
experimental  values  whereas  the  SA  solution  predicts  the 
breakdown  location  downstream  of  the  flight  test  and 
experimental  results.  Third,  the  SADES  solution  gives  a 
mean  value  of  vortex  breakdown  location  well  within  the 
flight  test  and  experimental  data.  It  should  also  be  noted 
that  the  computed  nondimensional  primary  frequency  of 
the  breakdown  oscillation  is  0.2  in  the  range  of 
frequencies  commonly  found  in  the  literature1 161  for  vortex 
breakdown.  This  inability  of  commonly  used  turbulence 
models  to  accurately  compute  a solution  with  breakdown 
is  well  documented  in  the  literature  and  is  due  to  the  large 
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amount  of  eddy-viscosity  these  models  put  into  the  core 
of  vortices.  Several  researchers  have  proposed  fixes  to 
these  turbulence  models  by  incorporating  some  form  of  a 
rotation  correction.  The  disadvantage  of  this  approach  is 
the  fact  the  simulation  will  still  be  operating  in  a RANS 
mode  and  compute  solutions  that  are  relatively  steady 
post-breakdown  as  opposed  to  an  LES  approach  that 
resolves  the  eddies  that  produce  the  unsteadiness.  It  is 
clear  in  Figures  5 and  6 that  the  SADES  method  does  not 
suffer  from  the  same  problem  as  the  RANS  methods  due 
to  the  fact  that  eddy  viscosity  is  computed  based  on  sub- 
grid scale  turbulence,  automatically  minimizing  the 
amount  of  spurious  eddy-viscosity  that  is  placed  in  the 
core  of  vortices. 

Figure  7 is  a well  known  plot  in  the  literature  of  the 
streamwise  location  of  the  LEX  vortex  breakdown  as  a 
function  of  angle-of-attack[l51.  The  current  solutions  fall 
in  the  range  of  flight  tests  and  experiments  plotted  at  30° 
angle-of-attack.  The  previous  comparisons  of  the  method 
with  the  flight  test  and  experimental  data  was  poor  due  to 
the  incorrect  flap  settings  and  diverter  slot  being 
uncovered.1'31 

This  section  presents  comparison  of  the  computed 
SADES  solutions  with  F- 1 8 HARV  flight  test  data  from 
NASA  Dryden.  The  HARV  was  instrumented  with  32 
kulite  pressure  sensors,  half  on  the  inboard  and  half  on  the 
outboard  sections  of  the  right  vertical  tail  (Figure  8).  The 
kulite  pressures  were  stored  every  30ms  as  a function  of 
time.  The  available  pressures  were  stored  relative  to  a 
reference  pressure  that  is  unfortunately  unknown.  The 
lack  of  known  reference  pressures  allowed  only  frequency 
comparisons  rather  than  frequency  and  amplitude 
comparisons  of  the  SADES  data  with  flight  test  data. 
Pressure  ports  of  Figure  8 circled  in  red  are  those  used  for 
comparison  with  the  SADES  simulations. 

The  flight  test  and  SADES  simulation  port  pressures 
were  analyzed  with  MATLAB’s  PSD  function.  Since  the 
flight  test  data  has  a different  time  step  and  period  of  time 
(40  sec),  the  power  resulting  from  a PSD  analysis  will  not 
be  a one  to  one  match  but  the  frequencies  and 
characteristic  shapes  of  the  PSD  should  match.  All  32 
pressure  ports  were  analyzed  but  only  a representative  set 
are  shown.  Figures  9a-b  depicts  the  comparison  of 
SADES  and  flight  test  data.  Figure  9a  shows  the  PSD 
data  for  flight  test  and  SADES  simulation  for  ports  17  and 
18  and  9b  shows  ports  25  and  26.  In  all  cases,  the 
frequency  content  shows  quite  good  comparison  between 
the  flight  test  and  SADES  simulations.  All  of  the  ports 
show  a wide  peak  amplitude  range  corresponding  to 
Strouhal  numbers  between  0.45  and  0.8  for  both  flight  test 
and  SADES  simulations.  This  frequency  range 
corresponds  to  pressure  sweeps  over  the  tail  surface 
observed  in  a movie  clip  of  the  SADES  simulation. 
Unfortunately,  the  published  first  bending  mode  is  at  a 
Strouhal  number  of  approximately  0.66  explaining  why 


the  tail  is  so  aeroelastically  active  at  this  flight  condition. 
Most  of  the  ports  also  show  matches  in  slopes  of  the  PSD 
for  the  Strouhal  range  of  1-10.  It  is  also  interesting  to 
note  that  when  the  flight  test  curves  for  each  port  lie  on 
top  of  each  other  this  is  true  for  the  SADES  solutions  as 
well  (Figures  9a  and  b),  and  when  the  flight  test  curves 
are  separated  they  are  separated  by  approximately  the 
same  amount  in  the  SADES  solutions  (Figure  9b).  A 
consistency  is  noted  in  the  level  of  power  between 
inboard  and  outboard  ports  for  both  flight  test  and 
SADES,  i.e.,  when  the  inboard  port  has  a1  higher  power 
for  flight  test  that  is  true  as  well  for  the  SADES 
simulation.  Finally,  when  the  curves  cross!  this  occurs  at 
approximately  the  same  frequency  for  flight  test  and 
SADES  (Figure  9b).  The  overall  comparison  of 
frequency  content  is  remarkably  good  for  the  SADES 
solutions,  demonstrating  the  utility  of  the  method  for  tail 
buffet  computations  at  flight  Reynolds  numbers. 

The  last  configuration  analyzed  was  in  direct 
response  to  a DoD  customer  request1  to  help  in 
understanding  the  nonlinear  aerodynamics  causing 
undesirable  flight  mechanics  of  a current  vehicle  in  a 
flight  test  program.  The  ARGUS  missile  experiences 
coning  after  release  making  the  trajectory  difficult  to 
predict.  The  simulations  compared  favorably  with 
experiments  and  a greater  understanding  of  the 
undesirable  phenomena  has  been  obtained.j  The  ARGUS 
missile  simulations  are  preliminary  to  six  degree  of 
freedom  simulations  of  the  free  flight  of  the  missile 
following  release. 

i 

3.4.  ARGUS.  i 

I 

The  ARGUS  program  is  intended  to!  meet  the  Air 
Force  requirement  of  detecting,  tracking,  identifying  and 
reporting  Time  Sensitive  Targets  in  near-real-time.  The 
first  generation  version  of  ARGUS  (Steel  Eagle  II)  had  an 
asymmetric  geometry  which  resulted  in  undesirable 
stability  characteristics.  A newly  designed  ARGUS 
projectile  was  recently  created  by  Textron,  Inc.  to  achieve 
improved  aerodynamic  characteristics.  To  complement 
the  development  effort,  CFD  analysis  of  the  new  design  is 
being  conducted  by  the  Academy’s  Department  of 
Aeronautics  to  determine  the  projectile’s  jlift,  drag,  and 
aerodynamic  moment  characteristics.  Problems  that  were 
encountered  with  the  previous  Steel  Eagle  II  design 
included  instability  in  the  sensor  air  body  when  deployed 
from  a carrier  aircraft  or  helicopter,  coning  instability 
during  free  flight,  non-zero  impact  angles,  and  high 
impact  velocities.  Most  of  the  problems  mentioned  above 
relate  to  the  coning  motion  of  the  vehicle.  While  it  was 
apparent  that  the  vehicle  displayed  coning  in  flight, 
neither  flight  testing  nor  wind  tunnel  testing  could  fully 
discover  why  the  motion  was  taking  place.  DES  of  the 
flow  field  for  the  ARGUS  geometry  was  conducted  at  M 
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= 0.5  and  a = 0°  and  vortex  shedding  was  detected  from 
the  terra  brakes  (Figure  12).  Researchers  at  USAFA 
suggested  that  the  vortex  shedding  could  be  alleviated  by 
drilling  holes  into  the  terra  brakes,  so  a DES  simulation  of 
that  configuration  was  also  conducted  at  the  same  flight 
conditions  (Figure  12).  While  the  flow  field  behind  the 
terra  brakes  with  holes  is  still  unsteady,  all  evidence  of 
vortex  shedding  is  gone — the  resulting  lateral  forces  and 
moments  have  shown  a corresponding  reduction  due  to 
the  addition  of  the  brakes,  which  should  alleviate  the 
coning  problem.  Initial  drop  tests  from  a helicopter  have 
verified  these  results,  and  the  ARGUS  geometry  will  be 
designed  with  holes  in  the  terra  brakes. 

3.  Conclusions 

The  proposed  method  of  solution  was  used  in  both 
basic  and  applied  simulations  in  the  Challenge  C92 
project  during  FY04.  The  basic  research  effort  made 
preliminary  progress  in  broadening  the  application  of  the 
method  to  laminar-turbulent  transition  and  embedded 
LES.  The  applications  made  great  strides  towards  aiding 
flight  test  of  full  aircraft  in  the  most  difficult  portions  of 
their  operational  envelopes.  The  F/A-18E  and  F/A-18C 
showed  excellent  comparison  with  experiment  and  flight 
tests  lending  credibility  to  the  method.  In  addition,  the 
method  was  used  to  aid  in  a current  DoD  flight  test 
program,  ARGUS,  by  helping  to  understand  complex 
nonlinear  aerodynamics  observed  in  flight  test  and 
experiment  resulting  from  massively  separated  flow.  The 
ARGUS  missile  system  simulations  resulted  in  a design 
change  of  the  vehicle. 
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Figure  1.  Contours  of  instantaneous  voracity  in  four 
planes  along  the  Aerospatiale-A  airfoil  at  an 
angle-of-attack  of  13.3  degrees  and  a Reynolds 
number  of  2*10® 
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Figure  2.  Normal  force  vs.  angle  of  attack  for  near  zero 
sideslip  for  the  F/A-18E  at  M=0.9 


Figure  3.  Rolling  moment  (left)  and  yawjng  moment 
(right)  vs.  bank  angle  for  the  F/A-18E  at  7°  pitch  angle 


Figure  4.  DES  forced  oscillations  simulations.  Left 
pane,  instantaneous  flow  visualization  [at  6°  pitch 
angle  - contours  of  pressure  on  surface,  and 
isosurface  of  zero  streamwise  velocity  (grey).  Right 
pane,  rolling  moment  vs.  roll  rate  phase  averaged  over 
5 cycles. 
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a)  SST  b)  SA  c)  SADES 

Figure  5.  Isometric  views  of  the  F/A-18C  at  a = 30°, 
Rec=  13*10®,  leading  edge  flaps  set  to  -33°,  trailing 
edge  flaps  set  to  0°,  with  no  diverter  slot  present:  a) 
SST  turbulence  model,  b)  SA  turbulence  model,  and  c) 
SADES  turbulence  model. 


Figure  6.  Time  histories  of  the  streamwise  coordinate 
of  vortex  breakdown  referenced  to  the  vehicles  nose 
and  scaled  by  the  length  for  the  SST,  SA,  and  SADES 
methods 
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Figure  7.  Streamwise  LEX  vortex  breakdown  position 
as  a function  of  angle-of-attack,  extracted  from 
Reference  47  SADES  mean  vortex  breakdown  position 
in  red. 


Figure  8.  Placement  of  the  F-18  HARV  Kulite  pressure 
sensors  on  the  right  vertical  tail.  Odd  port  numbers 
are  on  the  inboard  section  of  the  tail  and  even  are  on 
the  outboard  section.  Red  circles  around  ports 
indicates  those  used  in  comparing  flight  test  to 
SADES. 


from  Flight  Test  and  DES  Prediction  for  a)  Ports  17  & 
18,  and  b)  Ports  25  & 26 


Figure  10.  Time-averaged  induced  power  as  a function 
of  the  vertical  velocity  for  a single  rotor  with/without  a 
symmetric  boundary  condition,  and  a V-22  half  aircraft 
and  full  V-22  aircraft 
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Figure  11.  a)  Instantaneous  cross-plane  of  vorticity 
contours  for  the  V-22  in  hover  and  b)  Instantaneous 
cross-plane  of  vorticity  contours  for  the  V-22  in  a 
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Figure  12.  Iso-surfaces  of  vorticity  for  the  ARGUS 
simulations  for  terra  brakes  with  and  without  holes  at 
a Mach  number  of  0.5  and  an  angle-ofjattack  of  0 
degrees 


i 


I 

i 


i 


I 

I 


111 


